Reaching optimally oriented molecular states by laser kicks 
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We present a strategy for post-pulse molecular orientation aiming both at efficiency and maximal 
duration within a rotational period. We first identify the optimally oriented states which fulfill both 
requirements. We show that a sequence of half-cycle pulses of moderate intensity can be devised for 
reaching these target states. 

PACS numbers: 33.80.-b, 32.80.Lg, 42.50.Hz 



Molecular orientation plays a crucial role in a wide 
variety of applications extending from chemical reaction 
dynamicSjto surface processing, catalysis and nanoscale 
design 0, 0, 0, Q . Static electric field || and strong non- 
resonant long laser pulses 0, Q have been shown to yield 
adiabatic molecular orientation which disappears when 
the pulse is off. Noticeable orientation that persists after 
the end of the pulse (and even under thermal conditions) 
is of special importance for experiments requiring field- 
free transient orientation. It has recently been shown 
that very short pulses combining a frequency u) and its 
second harmonic 2tu excite a mixture of even and odd 
rotational levels and have the ability to produce such 
post-pulse orientation But even more decisive has 
been the suggestion to use half-cycle pulses (HCPs) , that 
through their highly asymmetrical shape induce a very 
sudden momentum transfer to the molecule which orients 
under such a kick after the field is off 0, ^| . Both the 
(lu + 2lo) and the kick mechanisms have received a con- 
firmation from optimal control schemes . The caveat 
is that the post-pulse orientation is maintained for only 
short times. Recently, the use of a train of kicks to in- 
crease the efficiency of the orientation has been suggested 
in optimal control strategies and applied to molec- 
ular alignment ^3 an d orientation of a 2D rotor [l3j . 
However, due to the strength of the kicks used, only the 
efficiency of the process has been optimized, its duration 
decreasing strongly. In the present letter, we propose a 
control strategy using specially designed series of kicks 



delivered by short HCPs, that allows to significantly en- 
hance the duration of the orientation, maintaining a high 
efficiency. Our construction is first based on the identifi- 
cation of target states which fulfill the previous require- 
ment. These states are characterized by the fact that 
they only involve a limited number of the lowest lying 
rotational levels and that they maximize the orientation 
efficiency within the corresponding restricted rotational 
spaces. At a second stage, we show that these selected 
states can be reached by a train of kicks, acting at appro- 
priately chosen times. The choice of the strength of the 
pulses (taken equal for simplicity), together with the to- 
tal number of kicks allow to approach these target states 
with good accuracy. 

The time evolution of the molecule (described in a 
3D rigid rotor approximation) interacting with a lin- 
early polarized field is governed by the time-dependent 
Schrodinger equation (in atomic units) 



[BL 2 - fj, F(t) cos 8]ip(0,<j),t), (1) 
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where L is the angular momentum operator, B the ro- 
tational constant, /zo the permanent dipole moment and 
F(t) the field amplitude. 8 denotes the polar angle be- 
tween the molecular axis and the polarization direction 
of the applied field. The motion related to the azimuthal 
angle <j> can be separated due to cylindrical symmetry. 
From now on, we assume a sudden approximation due to 
the short durations r of the HCPs, as compared to the 
molecular rotational period T rot = tt/B. For relatively 
low I (where I labels the quantum eigenstates of L 2 ), this 
amounts to the definition of a dimensionless, small per- 
turbative parameter e = tB. This definition, together 
with a rescaling of time s = t/r (such that s e [0,1] 
during the pulse) leads to an equation suitable for the 
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application of time-dependent unitary perturbation the- 

ory mm 
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[el? - E(s) cos 6]ip(d,(t>,s), (2) 



where E(s) = [jlqtF(ts). At lowest order in e, the emerg- 
ing dynamical picture is the following m^ 13 : an in- 
dividual HCP imparts a kick to the molecule described 
by an effective instantaneous evolution operator e lAcos0 j 
where A = E(s)ds is the total pulse area. Between 
two kicks, the molecule evolves under the effect of its 
field- free rotation e~ lsL s . 

The goal of the field driven molecular orientation is to 
maximize (or to minimize, depending on the choice of 
the orientation) for the longest time duration, the expec- 
tation value (cos0)(s) = (ip(6, <j>, s)\ cos 6\tp(6, 0, s)) after 
the pulse is over. An understanding of this process can 
be obtained by analysing the molecular dynamics in a fi- 
nite subspace Hm^ generated by the first (N + 1) eigen- 
states of L 2 , i.e. \l, m) (I — \m\, \m\ + 1, ...N + \m\) for 
a molecule initially in the state |/o > m,m). The justi- 
fication of such a reduction is in relation with the finite 
amount of energy that a finite number of HCPs of a given 
area A can transfer to the molecule. The mathematical 
advantage it offers is the consideration of an operator 



C^ = pW coa9P W t 



(3) 



{Pm^ being the projector on the subspace TLm ^) which, 
as opposed to cos has a discrete spectrum. It turns 
out that the state Iy^' 1 ), which maximises the orien- 



tation in the subspace Hm , is the eigenstate of C, 1 , 
with the highest eigenvalue. Using the approximation 
{l,m\cosd\l ±1, to) ~ 1/2 valid for I 3> to, straightfor- 
ward algebra leads to 
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(4) 

Since the quantum number to, related with the azimuthal 
angle is conserved, we will not write it explicitely, unless 
necessary. The maximal orientation in this subspace is 
found to be 
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(5) 



For a temperature T = K (m = 0), Fig. ^ gath- 
ers two informations relevant for the characterization of 
orientation as a function of N; namely (x^l cos^jx^') 
which is the maximum efficiency [approximatively given 
by Eq. J^J] that can ideally be expected for a process 
that stays confined within the finite subspace liS N > , and 
Ai/T rot which measures the relative duration of the ori- 
entation over which (cos 9) remains larger than 0.5 during 
the field-free evolution of \x" N '}- The results, expressed 
as a fraction of the rotational period T rot , are molecule 




FIG. 1: Maximal orientation efficiency (crosses) and associ- 
ated duration (open circles) as a function of N, where iV + 1 is 
the dimension of the rotationally excited subspace Ti.^^ (see 
text). The solid and dashed lines are just to guide the lecture. 



independent. From Fig. ^ we observe that, in order to 
keep a duration of the order of 1 /w of the rotational pe- 
riod (which may amount to durations exceeding 10 ps, 
for heavy diatomics like e.g. Nal), N has to be limited 
to 5 or 6, which seems rather limiting. But this turns 
out to be sufficient for very efficient orientation. N = 4 
already allows an orientation efficiency larger than 0.91. 

Two basic questions are in order : which set ot pa- 
rameters and which number of kicks have to be chosen 
to (approximately) remain in the subspace 1iS N \ and 
which strategies have to be followed to reach the maxi- 
mum possible efficiency within this subspace. The first 
question is in relation with the kick momentum transfer 
operator e lA cos 6 , which is the only evolution operator (as 
opposed to the free evolution) that rotationally excites 
the system and forces it to expand on a larger subspace. 
We can estimate the loss outside assuming a pre- 

liminar convergence to |x ) ( as shown by Fig. 0) by 
looking for the smallness of the norm 



— p( N ) e iAcose p( 



which, for small A, amounts to 



N' 



AN)\ 2 _ 



(AnY 



2(A + 2) 



(6) 



(T) 



This allows us to establish a relation between A and N 
for a given loss. N ~ 4 is found compatible with an 
i] ~ 0.02 (not more than 2% of the rotational population 
leaving the subspace 7i' 4 ') as far as A does not exceed 1 . 

The second question can be answered by adapting the 
strategy suggested for the orientation of a 2D rotor in 
Ref. [Fj> which consists in applying laser pulses each 
time (cos 0) reaches its maximum. The following argu- 
ment shows that, if the dynamics stays within the sub- 
space H^ N ', such a strategy precisely converges to an 
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optimal state |x^)- This is done by approximating the 
operators cos# and e lAcos6 by and e tAciN) respec- 

tively. The interaction with a sudden HCP only alters 
the slope of (C^ N ')(s) and not its value as is clear from 
the following relation 

L-iAC^ C {N) e iAC^ = / C (N)\ (8) 

Moreover, if a sudden pulse is applied at a time Si when 
<C (JV) }0) reaches its maximum d = (C^ N ^)(si), the 
slope undergoes a change from zero to a finite value 

d-(CW)\ Si+0 = i(e-* Ac(N) [eI 2 ,CW]e !Ac(, ")^0. 

(9) 

As (C^) is a periodic, continuously differentiable func- 
tion, it will reach within the rotational period, a max- 
imum value larger than the one obtained prior to the 
application of the pulse. Iterating the strategy, we get 
an increasing but bounded and therefore convergent se- 
quence of Ci's. Its limit is a fixed point C; = Cj+i, corre- 
sponding to the eigenvectors of the impulsive propagator 
e iAC {N) w i]jch are also the ones of C^ N K Indeed, at a 
fixed point, the slopes before and after the interaction 
with the last pulse have to be zero. 

Figure|2]gives two different views of the orientation dy- 
namics under the effect of a train of HCPs, separated by 
time delays corresponding to the above discussed strat- 
egy of maxima, with identical durations e = 0.01 and 
pulse areas A = 1 (that is about 0.3 ps and a field am- 
plitude of 1.5 ■ 10 5 V cm" 1 for LiCl 9]), which leads to 
a dynamics that remains within the subspace Tt^ for 
the considered numbers of kicks. From panel (a) it is 
interesting to note that a single kick produces an orien- 
tation of about 0.5, whereas the appropriate application 
of 15 kicks increases this efficiency up to 0.89, which is 
almost the optimal limit as found from Fig. ^ O n the 
other hand, the comparison of the average (cos 9) calcu- 
lated with the exact wave function \ip(s)) and with the 
one propagated in the subspace , are close enough to 
support the claim that the rotational dynamics actually 
resides within . Panel (b) shows the way the wave 
function \fp(s)) gets close to the optimally oriented state 

Ixl) 4 ''), showing thus the succesfull outcome of the pro- 
cess. The difference between the two dynamics can be 
also estimated by Eq. Q). Here again, the close conver- 
gence of the two calculations shows a coherent choice of 
A, N and the number of kicks for appropriately describ- 
ing the dynamics. The post-pulse dynamics, which is our 
main concern, is displayed on Fig. [3] [panel (a)]. A result 
expected from the previous analysis, but particularly re- 
markable with respect to previous proposals, is obtained 
with an efficiency of about 0.89 and a duration of the 
order of 2 /io of the rotational period (that is about 2 ps 
for a light molecule like LiCl and 20 ps for a heavy one, 
like Nal). 




rot 

FIG. 2: Orientation dynamics during the train of HCPs at 
T = K: panel (a) for {ip(s)\ cos 9\tj)(s)) and panel (b) for 

| (Xq 4 " 1 ^(s)) | 2 - The solid line corresponds to \ip(s)) calculated 
exactly and the dashed line to the wave function propagated 
in the subspace . The train of HCPs is displayed on panel 
(b) and the optimal orientation is indicated by an horizontal 
line on panel (a). 

Orientation is subject to a drastic decrease with tem- 
perature [HI El- This is basically due to the fact that, 
with non-zero temperature, the initial state is a superpo- 
sition of a statistical ensemble of rotational states with 
to 7^ 0, which tends to misalign the molecule. The effi- 
ciency of the orientation is characterized by an additional 
average of (cos#) over the density operators p m {s) 

((cos0))(s) = 1*M») cos 0]. (10) 

We recall that p m (s) evolves according to the von Neu- 
mann equation 

^Pm(s) = i[p m {s),eL 2 - E(s) cos0], (11) 
with as initial condition 

Pm(0) = | E \l,m)e- Bl{ - l+1)/kT {lM, (12) 

l>\m\ 

where Z = E m ezEi>\ m \^ mi+1)/kT is the partition 
function and k the Boltzmann constant. Following our 
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FIG. 3: Post-pulse orientation dynamics of the molecule LiCl 
after interaction with a train of HCPs (the time t — corre- 
sponds to the first kick) in the case A = 1, T = K [panel 
(a)] and A — 2, T — 5 K [panel (b)]. Solid and dashed lines 
correspond respectively to the averages calculated with the 
exact wave function and the optimal state [|Xo > on panel 
(a) and p' 7 ' on panel (b)]. 

previous analysis, we are looking for the optimal den- 
sity operator which maximizes {{cos 9}} in the subspace 
H^ N \ which is given by = J2 m P^ where pm = 



\Xm )Tr[pm }{xm I with the constraint related to the 
conservation of m, expressed as Tr{pffl] to be kept con- 
stant 

l=\m\+N 

Tt[pW] = J_ £ e -«<H-i)/*r (13) 

|m| 

The effect of a sudden pulse on p is e~ lA cos 6 pe lA cos 8 , and 
because its optimal value corresponds to a fixed point of 
the {{C^ N ^))(Si) sequence, it is precisely the application 
of a train of HCPs with individual pulses at times Si 
where ((cos#))(sj) reaches its maximum, that converges 
to the best possible orientation within this model. The 
resulting dynamics is plotted in Fig. [panel (b)] , with 
a maximum efficiency of about 0.75 and a duration of 
about 1 /2o of the rotational period. To our knowledge, 
this is the largest duration and efficiency achieved up to 
date for a thermal ensemble. 

In conclusion, we have presented tools for controlling 
molecular orientation dynamics using a train of HCPs, 
achieving both efficiency and duration of the post-pulse 
orientation. Moreover, this scheme can be expected to 
be transposable to a generic system, with free periodic 
dynamics governed by an Hamiltonian Ho , and for which 
we are aiming to optimally control an observable O (i.e. 
maximize or minimize the average (0}(t) of an upper 
or lower bounded operator O which does not commute 
with Hq ). This could be done through a device that 
perturbs the system according to a unitary operator U, 
which commutes with O, such that its application does 
not alter (O) = (JJ~ l OU) on one hand, and the optimal 
target state is an eigenfunction of both O and U on the 
other hand. This optimum corresponds to a fixed point 
of the sequence Oi = (0)(ti) where ti are the times when 
(O) (t) reaches its maximum (or minimum) under the free 
evolution. In particular, this scheme provides a compre- 
hension of previous works on alignment and orientation 
of 2D and 3D rotors [HEi. 



[1] P. R. Brooks, Science 193, 11 (1976). 
[2] F. J. Aoiz, Chem. Phys. Lett. 289, 132 (1998). 
[3] T. Seideman, Phys. Rev. A 56, R17 (1997). 
[4] H. Stapelfeldt and T. Seideman, Rev. Mod. Phys. 75, 
543 (2003). 

[5] L. Cai, J. Marango and B. Friedrich, Phys. Rev. Lett. 

86, 775 (2001). 
[6] M. J. J. Vrakking and S. Stolte, Chem. Phys. Lett. 271, 

209 (1997). 

[7] S. Guerin, L. P. Yatsenko, H. R. Jauslin, O. Faucher and 
B. Lavorel, Phys. Rev. Lett. 88, 233601 (2002). 

[8] C. M. Dion, A. D. Bandrauk, O. Atabek, A. Keller, H. 
Umeda and Y. Fujimura, Chem. Phys. Lett. 302, 215 
(1999). 

[9] C. M. Dion, A. Keller and O. Atabek, Eur. Phys. J. D. 



14, 249 (2001). 

[10] M. Machholm and N. E. Henriksen, Phys. Rev. Lett. 87, 
193001 (2001). 

[11] C. M. Dion, A. Ben Haj Yedder, E. Cances, A. Keller, C. 

L. Bris and O. Atabek, Phys. Rev. A 65, 063408 (2002). 
[12] M. Leibscher, I. Sh. Averbukh and H. Rabitz, Phys. Rev. 

Lett. 90, 213001 (2003). 
[13] I. Sh. Averbukh and R. Arvieu, Phys. Rev. Lett. 87, 

163601 (2001). 

[14] D. Daems, A. Keller, S. Guerin, H. R. Jauslin and O. 

Atabek, Phys. Rev. A 67, 052505 (2003). 
[15] D. Sugny, A. Keller, O. Atabek, D. Daems, S. Guerin 

and H. R. Jauslin, Phys. Rev. A submitted (2003). 
[16] N. E. Henriksen, Chem. Phys. Lett. 312, 196 (1999). 



